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Abstract. 

Self-consistent modelling based on local spin-density formalism is employed to 
calculate conductance of quantum point contacts at finite temperatures. The total 
electrostatic potential exhibits spin-dependent splitting, which persists at temperatures 
up to 0.5 K and gives rise to the anomalies at 0.4 and 0.7 of the conductance quantum 
(2e 2 /h) occurring simultaneously in the absence of external magnetic fields. 



PACS numbers: 73.61.-r, 71.15.Mb. 71. 70. Cm 



| Present address: Department of Physics, University of Queensland, Brisbane, Australia, 
stace@physics.uq.cdu.au. 



0.4 and 0.7 conductance anomalies in quantum point contacts 



2 



1. Introduction 

The quantisation of quantum wire conductance in multiples of Go = 2e 2 /h has long 
since been observed p~j [2] and understood [3]. The extreme sensitivity of short 
quantum wires formed in the constriction of quantum point contacts (QPCs) to the 
local electrostatic environment is demonstrated in their use as charge detectors capable 
of monitoring single-electron tunnelling in real time [HE]. In addition, the first 
conductance plateau exhibits a series of non- integer features, most notably at 0.7Go [6], 
observed in various semiconductor quantum wires under a range of experimental 
conditions [7j. It is largely agreed that the non-integer conductance features are due 
to electron-electron interactions of the confined 2DEG, consistent with the fact that 
the 0.7 structure evolves into a spin-split subband at 0.5G with an in-plane magnetic 
field [6], but a comprehensive understanding of the effect is yet to emerge. Proposals 
based on microscopic many-body theory ascribe the observed conductance behaviour 
to the spontaneous spin polarisation [El [9], dO, EEH EE] , and Kondo-like correlated spin 
state [TBI [H] in the QPC, whereas phenomenological models assume the pinning of one 
of the spin subbands to the chemical potential at the opening of a new conductance 
channel pH [HI H7] . 

In this paper, we use the local spin-density approximation (LSDA) model developed 
by Berggren et al. [UJ . The model is modified to include the case of finite temperatures, 
by replacing the zero-temperature Fermi distribution with the finite temperature version. 
We use the LSDA results to calculate the conductance using a rigorous form of the 
transmission coefficient obtained from the self-consistent electrostatic potential. We 
find that the transverse modes split into spin-dependent subbands, both above and 
below the chemical potential. This results in conductance features at OAGq and 0.7Gq, 
consistent with the latest experimental results of Crook et al. [18]. In our model, the 
features persist at temperatures up to 0.5 K, and agrees qualitatively with recent results 

2. Methods 

2.1. Self- consistent electrostatic potential 

Our simulated QPC device is a modulation-doped heterostructure, which consists of an 
AlGaAs spacer, an AlGaAs donor layer (uniformly doped at = 6 x 10 17 cm -3 ), and a 
GaAs cap grown epitaxially on a GaAs substrate. A plane metallic gate is deposited on 
top of the cap; two 100 x 400 nm 2 rectangular openings connected by a narrow channel 
200 nm long and 10 nm wide are then etched through the gate (see Figure (IK). Applying 
a negative voltage between the gate and the substrate creates electrostatic confinement 
of the 2DEG at the heterojunction 70 nm below the gate plane thus forming a QPC. By 
ensuring the constriction region is far from the boundaries, we may neglect boundary 
effects, as discussed in [9]. 

To obtain the total spin- dependent electrostatic potentials for the confined 2DEG, 
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we employ an LSDA theoretical framework described in detail elsewhere [TT j [51 [T51 12"0] . 
The robustness of this computational method, at zero temperature, has been established 
in these studies. The model takes into account electron exchange-correlation effects, 
whereby the exchange term U" oc —\fn? is a function of the 2D electron density, and 
the correlation term is taken in the parameterized form of Tanatar and Ceperley [21] 
as an interpolation between the fully spin-polarised and non-polarised cases [22] . Single 
particle LSDA equations for spin up and spin down electrons are solved self-consistently 
via an iterative procedure with the input values for the potentials derived from the 
semiclassical Thomas-Fermi approximation. For our system, we have a total of around 
300 electrons in the region being modelled (yielding an electron density of about 
8 x 10 10 cm -2 ), of which around 10 are in the constriction. 

Finite temperature is included by calculating the electron density at each iteration 

as 



where Ef and ipf(x,y) are the eigenenergies and eigenvectors of the LSDA equations, 
respectively, and f(E[, fx) = (1 + e( E i~ /J, >' kT )~ 1 is the Fermi-Dirac distribution function 
with the chemical potential set at fx = 0. 

2.2. Conductance 

Having computed the effective single particle potential, we can now calculate the 
conductance in the linear regime according to 



where Ty s (E) is the spin and energy dependent transmission coefficient, and f(E) 
is the Fermi-Dirac occupation. The dependence on an external gate voltage, V g , is 
included explicitly here, but will be dropped subsequently. To compute the transmission 
coefficients we make a ID approximation for the longitudinal potential: we consider a 
particle scattering from a ID potential given by the lowest lateral subband energy, 
determined from the self-consistent potential. The details of the calculation are 
described in the Appendix. 

3. Results 

Figure [1] shows a typical self-consistent total electrostatic potential and its longitudinal 
and transverse cross-sections at the saddle point. Spin-dependent splitting of the 
potentials occurs in the absence of an external magnetic field, driven mostly by the 
electron exchange term [8j [TTJ [9j [23]. In the following, we consider only the central 
300 x 400 nm 2 region of the device where the potentials are not affected by the finite 
size of the lithographic gate openings on either side of the channel. 




(1) 




(2) 
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Figure 1. (a) Top view of the QPC device showing lithographically defined openings 
surrounded by the metallic gate. Dashed line outlines the region of consideration, (b) 
Total electrostatic potential for spin up electrons at T = 0.1 K and V g = —0.5007 V. 
(c,d) Saddle-point longitudinal and transverse, respectively, cross-sections of the spin 
up and spin down potentials show spin gaps in the absence of the magnetic field. 




Figure 2. Difference between spin up and spin down potentials along (a) the 
longitudinal centreline, (b) in the lateral direction through the middle of the saddle, 
for different gate voltages at T = 0.1 K. The colour bars show the scale in meV. 



0.4 and 0.7 conductance anomalies in quantum point contacts 



5 



1.0 - 




0.8 - 



I 06 - 



£ 0.4 - 



CD 



0.2 - 



0.0 - 



-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
E-\i [meV] 



Figure 3. Transmission coefficient versus incoming energy at T = 0.1 K, at different 
gate voltages. The heavier lines are spin down and lighter lines are spin up. Energy 
is defined with respect to the chemical potential. Note that at V g = —0.4990^ the up 
and down traces are almost the same. 



Figure [2] shows the difference between the spin up and spin down potentials along 
the x and y directions, passing through the middle of the saddle. There are clearly 
two regions in which the potentials are split: in the region around V g = —0.504 V 
and V g = —0.501 V. This splitting is the origin of the features that we find in the 
conductance. Figure [2b shows the splitting in the lateral direction, demonstrating that 
the splitting is localised around the saddle-point. The constriction thus forms a quasi- 
1D channel. We note the applied gate voltages here differ in scale from experiments due 
to different device geometries. 

As described above and in the Appendix, we compute the transmission, T(E) as a 
function of the incident energy, which is used in Eq. (J2j) to calculate the conductance. 
Some representative examples are shown in Figure [3] for two different gate voltages. 
At V g = —0.501 V where the splitting is significant, the transmission coefficients are 
very different for the two spin configurations. At V g = —0.4990 V the splitting is very 
small, so the transmission coefficients are very similar. The width of the step in the 
transmission is a measure of the flatness of the top of the tunnel barrier; compare with 
Figure db, where the barrier vanishes over an energy scale of about 2 meV. Since the 
temperature for this example is 0.1K, the derivative of the Fermi-Dirac distribution in 
Eq. (j2j) is only ~ 0.01 meV wide. Over this energy scale, the transmission coefficients 
are essentially constant (see Figure [3]), so to a good approximation 



Substituting the numerical values from Figure [3] for T*(jj,) and T^(fi) at V g = —0.501 
gives G « 0.7G . 

The variation of the transmission probability with V g can be most readily 
understood by considering the dependence of the subband energies on V g , shown in fig.H] 
(right). The two splittings of the self-consistent potentials give rise to commensurate 



G = G [T T ( /U )+T i (/i)]/2. 



(3) 
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Figure 4. Subband energy of the lowest transverse mode for different gate voltages. 
The energies are defined with respect to the chemical potential. Lower curves are spin 
up subbands, and upper curves are spin down subbands. Note that there is almost no 
splitting at IK. 

splittings in the subband energies. The splitting around V g = —0.504 V does not 
track the chemical potential, but instead, the spin up subband dips below the chemical 
potential. This increases the transmission coefficient for the spin up component, giving 
rise to a peak in the conductance. Similarly, the spin down subband rises above the 
chemical potential around V g = —0.501 V, which suppresses the transmission coefficient 
for the spin down component, leading to a dip in the conductance. 

The two splittings give rise to two features in the conductance, shown in Figure at 
V g = -0.504 V and V g = -0.501 V respectively. Notably, the feature at V g = -0.501 V 
is strongly reminiscent of the 0.7Go structure seen in numerous experiments. There also 
a feature at 0.4Go seen simultaneously with the 0.7Go in the absence of external magnetic 
fields. Both features vanish with increasing temperature, the 0.4Go feature being slightly 
more robust than the 0.7Go- We are not sure of the origins of the small oscillations seen 
around the features at the lowest temperatures, but the general behaviour is robust |24j . 

4. Discussion 

The results presented above establish the relative insensitivity of the 0.7 feature to 
temperature, since it persists up to 0.5 K, for the model we considered. The reasons for 
the feature are largely compatible with phenomenological models pJ2 HEl EG] , with a few 
significant differences. In phenomenological model of [17] , the subband energy of one 
spin configuration dips below the chemical potential, whilst the other closely tracks the 
chemical potential over some finite range of V g . In contrast, our results indicate that the 
spin subbands do not split as they cross the chemical potential, but rather after they 
have fallen below it. After splitting, one subband (labelled "Down" in Figure Hj) rises 
through the chemical potential, thereby suppressing its transmission coefficient around 
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Figure 5. Conductance versus gate voltage, for different temperatures, 0.01 K, 0.05 
K, 0.1 K, 0.5 K and 1 K. The traces are shifted in voltage for clarity. 

V g = -0.501V. 

In contrast to many other theoretical and phenomenological investigations of 
this subject, we have computed the transmission coefficient from the self-consistent 
scattering potential. We found that it varies slowly on the scale of any temperature we 
consider, rather than behaving as a step function at the chemical potential. Furthermore, 
we are interested in the transmission just below the first plateau, which is when the WKB 
method breaks down. Thus, our treatment makes a fair estimate of T(E). 

Whilst the feature at OAGq was not entirely unexpected, having been discussed 
in [11], it is not commonly observed in experiment. Recent observations show a 
structure at around 0.4 to 0.5G7o [IB], which is ascribed to the mirror symmetry of 
the channel, consistent with the calculation presented here. This suggests a future 
avenue of exploration on the role of asymmetry in the channel. This is beyond the 
scope of the present work. Nevertheless, the 0AG$ feature arises in our calculation 
for the same reason as the O.IGq feature: the spin subbands split above the chemical 
potential, causing one of them to approach the chemical potential over a narrow range 
oiV g . 

In conclusion, we have computed the self-consistent scattering potential for a QPC 
at different temperatures, from which we calculated the transmission coefficient and 
conductance. We found features below the first conductance plateau consistent with the 
experimentally observed 0.7Go structure, as well as another at 0.467 . The features are 
robust against temperature, as seen in experiment, but the high temperature behaviour 
that we find is not in agreement with observation. 

Appendix: Scattering Calculation 

Consider a ID potential V(x) localised in the region \x\ < X. We compute the 
tunnelling amplitude for a particle entering the region from the left with energy E. 
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The wavefunction takes the form 

Ae ikx + Be~ ikx for x < -X 
iPe{%) = { <j){x) for \x\ < X (4) 

Ce lkx for x > X 

< 

where E = %k 2 E /2m, and both ip(x) and ip'(x) are continuous everywhere. The 
tunnelling probability is then given by T(E) = \C/A\ 2 . 4>{x) satisfies the ID Schrodinger 
equation in the scattering region, 

- h 2 /2m(f)"(x) + V{x)<j){x) =E<j>{x). (5) 

To satisfy the matching conditions at x = ±X, we first solve eq. ([5]) twice, using two 
linearly independent sets of boundary conditions, such as {(f>i(— X) = l,<f>i(X) = 0} 
and {02(— X) = 0, 4>2(X) = 1}, leading to two linearly independent solutions (pi and (j>2- 
These are then used to construct (p(x) = a<pi(x) + Pfoix), and the parameters B,C,a 
and (3 are solved to satisfy the four matching conditions. With the boundary conditions 
above, and for a symmetric scattering potential (so that faix) = <j>i(—x)), we find 

C 2te- 2ikX k(f)' 1 {X) 



A (fc-i#(-X)) 2 + &(X) 2 " 
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